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We determine the flow structure in an axisymmetric diffuser or expansion region connecting 
two cylindrical pipes when the inlet flow is a solid body rotation with a uniform axial flow of 
speeds O and U, respectively. A quasi-cylindrical approximation is made in order to solve 
the steady Euler equation, mainly the Bragg-Hawthorne equation. As in our previous work 
on the cylindrical region downstream [R Gonzalez et al., Phys. Fluids 20, 24106 (2008); 
R. Gonzalez et al., Phys. Fluids 22, 74102 (2010), R Gonzalez et al., J. Phys.: Conf. 
Ser. 296, 012024 (2011)], the steady flow in the transition region shows a Beltrami flow 
structure. The Beltrami flow is defined as a field vs that satisfies ujb = V x v_b = 7v_g, 
with 7 = constant. We say that the flow has a Beltrami flow structure when it can be put 
in the form v = Ue z +£lree + vb, being U and Q constants, i.e it is the superposition of a 
solid body rotation and translation with a Beltrami one. Therefore, those findings about 
flow stability hold. The quasi-cylindrical solutions do not branch off and the results do not 
depend on the chosen transition profile in view of the boundary conditions considered. By 
comparing this with our earliest work, we relate the critical Rossby number $ cs (stagnation) 
to the corresponding one at the fold § c f [J. D. Buntine et al., Proc. R. Soc. Lond. A 449, 
139 (1995)]. 
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I. Introduction 

We have recently conducted studies on the for- 
mation of Kelvin waves and some of their fea- 
tures when an axisymmetric Rankine flow expe- 
riences a soft expansion between two cylindrical 
pipes [TJ[2]. One of the significant characteristics 
of this phenomenon is that the downstream flow 
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shows a Rankine flow superposing a Beltrami flow 
(Beltrami flow structure [4 )). Yet, upstream and 
downstream cylindrical geometries were considered 
without taking into account the flow in the expan- 
sion. This work considered that the base upstream 
flow, formed by a vortex core surrounded by a po- 
tential flow, would have the same Beltrami struc- 
ture at the expansion and downstream. Neverthe- 
less, the flow at the expansion was not determined. 
However, it has been seen that this flow is only 
possible when no reversed flow is present and if its 
parameters do not take the values where a vortex 
breakdown appears [BHE]- The starting point in 
the study of the expansion flow is an axysimmetric 
steady state resulting from the Bragg-Hawthorne 
equation [71I5HTT] for both the vortex breakdown 
and the formation of waves. Therefore, the solu- 
tion behavior, whether it branches off or shows a 
possible stagnation point on the axis, will be deter- 



040002-1 



Papers in Physics, vol. 4, art. 040002 (2012) / R Gonzalez et al. 



minant to delimit both phenomena. 

Our previous research focused on the formation 
of Kelvin waves with a Beltrami flow structure 
downstream [TH3], when the upstream flow was a 
Rankine one. This present investigation considers 
only a solid body rotation flow with uniform axial 
flow at the inlet. As a first step in the study of the 
flow at the expansion, we only study the rotational 
flow. However, comparisons with our previous work 
[T] will be drawn. 

The aim of this present work is to obtain the 
steady flow structure at the expansion, consider- 
ing a quasi-cylindrical approximation when the in- 
let flow is a solid body rotation with uniform axial 
flow of speeds £7 and U, respectively. If a is the ra- 
dius of the cylindrical region upstream, a relevant 
parameter is the Rossby number $ — -M-. Thus, 
we would like to determine how this flow depends 
on the Rossby number, on the geometrical param- 
eters of the expansion and on the critical values of 
the parameters. We focus on finding the param- 
eter values for which a stagnation point emerges 
on the axis, or for which the solution of the Bragg- 
Hawthorne equation branches off. We take them as 
the conditions for the vortex breakdown to develop. 

First, this paper presents the inlet flow and the 
corresponding Bragg-Hawthornc equation written 
for the transition together with the boundary con- 
ditions in section II. Second, it works on the 
quasi-cylindrical approximation for the Bragg- 
Hawthorne equation and its solution is developed 
in section III. Third, results and discussions are of- 
fered in section IV together with a comparison with 
our previous work [T|. Finally, conclusions are pre- 
sented in section V. 



Euler equation which can be written as the Bragg- 
Hawthorne equation [TU] 



d 2 ^ d (Idij) 
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dz 2 dr \r dr 
where ip is the defined stream function 



I dip 

r dz ' 



I dip 
r dr ' 



(3) 



and H(ip), C(ip) are the total head and the circula- 
tion, respectively 



H{i>) = \{vl + v 2 + v 2 z ) + ^, C(V) = rvg. (4) 
2 p 

To solve Eq. @, the boundary conditions must 
be established. These consist of giving the inlet 
flow, of being both the centerline and the bound- 
ary wall, streamlines, and of being the axial veloc- 
ity positive (v z > 0). For the upstream flow, the 
stream function is ip — ^Ur 2 , and H(ip), C(tp) are 
given by 



W) = ^u 2 + n^, c(v>)=7^ (5) 

7 = 2M- being the eigenvalue of the flow with Bel- 
trami structure 3 . Thus, by considering the inlet 
flow, Eqs. ([5]) are valid for the whole region. The 
second condition regarding the streamlines implies 
the following relations 



II. The Bragg— Hawthorne equation 

We assume an upstream flow in a pipe of radius 
a as an inlet flow in an axisymmetric expansion of 
length L connecting to another pipe with radius b, 
b > a. The inlet flow filling the pipe consists of a 
solid body rotation of speed £1 with a uniform axial 
flow of speed U: 



= Ue z 



Clreg, 



(1) 



U and £1 being constants. The equilibrium flow 
in the whole region is determined by the steady 



V>(r = 0,z) = 0, 
ip(r 



a(z),z) = -Ua 2 , 



< z< L 



(6) 



where r = cr(z) gives the axisymmetric profile of 
the pipe expansion. Deducing from Eq. ©, the 
boundary conditions are determined by the inlet 
flow. Additionally, curved profiles are considered, 

so 



-8~z^ Z 



L) = 0, < r < b. 



(7) 
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III. Quasi-cylindrical approximation ii- curved profile 



r) 2 / - 

If we consider that -^ — 0, the solutions to Eqs. 
([2]) and ([5]) for the cylindrical regions are given by 



i 



ip = -Ur 2 +ArJi[yr}, 



(8) 



where A is a constant. The quasi-cylindrical ap- 
proximation consists of taking the dependence of 
A{z) on z but with the condition -^ k, com- 
pared with the remaining terms of (J2J . The ampli- 
tude A(z) is then obtained by imposing the bound- 
ary conditions © which depend on the wall profile 
r = o-(z), giving 



A(z) 



1U i 

2" 



a 2 (z)) 



(9) 



cr(z)Ji[7cr(z)] ' 

By using the dimensionless quantities r = - 
z = -, v = tj- the stream function in the quasi 
cylindrical approximation can be written as 



< z < L. 



(15) 



The latter meets the boundary condition ([7]) as 
well. Therefore, Eqs. (J1HI15P together with the 
boundary conditions (|6l7p allow to determine the 
flow structure for both the conical and curved wall 
profile. 

IV. Results and discussion 

We note that the flow keeps a Beltrami flow struc- 
ture in the quasi-cylindrical approximation. Effec- 
tively, giving (UJJ|13J) 



V r (r,z) 



v Br {f,z) 



(16) 



1 2 

V'= -f 2 +A(z)fJi[-f], 



A(z) 



1 (1-^)) 
2a(z)J 1 [^Cz)Y 



(10) 



where $ = jy- is the Rossby number. Hence the 
velocity field becomes 



v r (f,z) = -A(z)Ji[-f\ (11) 

Mf,z) = if+^i^Ji^f] (12) 

v*{f,z) = l + ^A(z)Jo[^r], (13) 

where A (z) = dA(z)/dz. 

Finally, it is necessary to give the wall profile 
ct(z) to completely determine the flow. Two kinds 
of profiles were seen: 

i- conical profile 



<->■( I) = 1 + ( ^-=— I :. 



< z < L and r\ = — . 
a 



(14) 



ve(r,z) 



f + v B e(f,z) 



v z (r,z) = l + VBz(f,z), 



(17) 



(18) 



it is easy to see that under this approximation V x 
Vs(f, z) = §Vs(f, z) and so, the whole flow is the 
sum of a solid body rotation flow with a uniform 
axial flow plus a Beltrami flow, given the latter in 
a system with uniform translation velocity U = l.z 
and uniform rigid rotation velocity V = hfd. 

Given the flow field and its structure, the param- 
eters are considered by evaluating the behavior of 
v z (f, z~o) with z~q — L i.e., taken at outlet, and with 
L = 1. In order to do so, a wall profile is selected 
(fT4l or [j"5|) and three different values of the expan- 
sion parameter are taken, mainly r\\ = 1.1, 772 = 1-2 
and 773 = 1.3. 

The first step is to analyze the flow dependence 
on the Rossby number. In Fig. [TJ the contour 
flows corresponding to the conical and curved pro- 
files for rji = 1.1, $1 = 0.695 are shown. Graph- 
ics in Fig. [5] represent the same configuration but 
for 1} — 0.68 < i?i. The broken lines represent 
points for which v z = 0. Inflow and recirculation 
are present but it is not a real flow because the 
model fails when considering inflow. It can be seen 
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Figure 1: Contour flow in the transition region for conical and curved profiles for r\\ — 1.1, #i = 0.695. 
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Figure 2: Contour flow in the transition region for conical and curved profiles for r/i = 1.1, i?j = 0.68. 
The broken lines represent points with v z = 0. 



that for i?i = 0.695, v z — at the outlet, on the 
axis. For the Rossby numbers with $ > i? c , the 
azimuthal flow vorticity is negative {uj^ < 0), re- 
sulting in an increase in the axial velocity with the 
radius, and so having a minimum on the axis where 
the stagnation point appears [BJ. Therefore, the 
critical Rossby number can be defined i? c as the 
value where v z is zero at the outlet on the axis i.e., 
where the flow shows a stagnation point. This is 



the necessary condition to produce a vortex break- 
down [6j . We find the same critical Rossby number 
for both wall profiles and so we will not treat them 
separately from now on. The critical Rosssby val- 
ues for rj2 = 1.2 and 773 = 1.3 are $2 = 0.869 and 
$3 = 1.052, respectively. 

Given the previous analysis, the second step is to 
show the behavior of v z on the axis at the outlet as 
a function of 1? for each 7/ in order to study the ex- 
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istence of folds in the Rossby number-continuation 
parameter (equivalent to the swirl parameter in 
[51171111] ); indeed, we have seen that v z has the min- 
imum on the axis. Besides, when using Eq. (|13[) 
when r = 0, it is easy to see that v z decreases with 
z and so it reaches the minimum at the outlet being 
v z > 0. In Fig. |31 the radial dependence of v z is 
plotted at the outlet for 771,772,773 and its variation 
with 7? when it is slightly shifted from #1. In Fig. 
[4J it can be seen that the minimum of v z on the 
axis increases with 7? so there is no fold of v zm i n 
as defined by Buntinc and Saffman in a similar ap- 
proximation [5J. 



n\<r]2< m 




(a) 




(b) 



Figure 3: (a) v z at the outlet as a function of r 
for 771,772,773 and the corresponding critical Rossby 
numbers $i,i?2)$3' 3 ) ^z a t the outlet as a function 
of r for $1 and for values of t? slightly shifted from 
t?i . In each case, the minimum of v z is reached on 
the axis. 



The dependence of the results on L is analyzed. 
It can be seen that when z — L in Eqs. (|14p and 
([j~5j) . a{L) — 77 is obtained. By replacing this in 
Eq. ([TB")) for z = L and r = it gives 



v,min = 1 



%/l [§*?]' 



(19) 



and so t? c is obtained as a function of 77 by solving 
the last equation when v z min = 0, as shown in 
Fig. [5j This result seems to be surprising, but 
it is not so if it is considered as derived from the 
quasi-cylindrical approximation: the dependence of 
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Figure 4: v z at the outlet on the axis as a function 
of the Rossby number $ for 771 = 1.1, 772 = 1.2, 773 = 
1.3. Here i?i = 0.695, t9 2 = 0.869 and t? 3 = 1.052 
correspond to stagnation points. 



the flow on z is obtained through the boundary 
conditions expressed by Eq. §§§. At the same time, 
these boundary conditions depend on the inlet flow 
and on the parameter 77. This explains the fact 
that the same results, for both conical and curved 
profiles, have been obtained and that the condition 
given by Eq. (J7J at the outlet has not influenced 
them. 
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Figure 5: Critical Rossby number d c as a function 


of 77. 





Differences with Batchelor's seminal work should 
be marked 40]. Mainly, he works in cylindrical ge- 
ometry and does not consider the dependence of the 
flow on z . We introduce this z dependence through 
the quasi-cylindrical approximation. This, there- 
fore, allows us to find the structure of the flow in the 
transition together with the Rossby critical number 
defined by considering this structure and by show- 
ing that the minimum of v z is reached at the outlet 
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on the axis. Nevertheless, once the flow reaches the 
pipe downstream, the analysis coincides because, as 
shown, the problem depends on the inlet flow and 
on the parameter expansion rj. This allows us to 
consider the issue of the vortex core that we have 
not considered at the inlet flow. As we know the 
structure of the flow in the downstream cylindri- 
cal region PQ and by assuming a quasi-cylindrical 
approximation for the vortex core in the transition 



region, the minimum of v c 
axis is given by 



at the outlet on the 



ucorezmrn 



1 + 



(1~7? 2 ) 



(20) 



where i? = — , f] = * and £ and l are the dimension- 
less radius of the core downstream and upstream, 
respectively. We note that r) is the expansion pa- 
rameter of the core. Hence Eqs. p^|) and (|2"0"|) 
have the same structure. In the present work, we 
have not found any fold in the Rossby number- 
continuation parameter of v z , as found in our pre- 
vious work PQ where the fold was associated with a 
critical Rossby number called f? c / by Buntine and 
Saffman [5]. As we have already done, we define 
the Rossby critical number for which v corezmin = 
where there is a stagnation point, and we will call 
it -& cs . In 1 , for l = 0.272 and pipe expansion 
parameters ^1,772,773, we have found that d c f were 
0.35, 0.44 and 0.53, respectively, while the core ex- 
pansion parameters f\ were 1.25, 1.47 and 1.65, re- 
spectively. 

By replacing these values in Eq. (|2"U|| when 

Vcorezmin * s zero j we get the corresponding i) cs and 
then # cs for the vortex core. These are respectively 
0.26, 0.38 and 0.49. That is to say that in all the 
cases we have $ cs < i? c /. Therefore, at the fold 
v z > 0. This coincides with the results found by 
Buntine and Saffman [5] in their analysis using a 
three-parameter family inlet flow. 

V. Conclusions 

The main conclusions drawn from the previous sec- 
tions are: 

1. In the quasi-cylindrical approximation, the 
steady flow in the transition expansion region 
corresponding to a solid body rotation with 



uniform axial flow as inlet flow has the same 
Beltrami flow structure as in the pipe down- 
stream, which is compatible with the bound- 
ary conditions. Therefore, findings from our 
previous work on stability [IH3| can hold. 

2. For fixed values of ij and 1? > i? c , lu^, < and 
then v z in the transition region is an increasing 
function of r and a decreasing function of z 
reaching its the minimum on the axis at the 
outlet. 

3. For fixed values of 77, the minimum of v z on 
the axis is an increasing function of d (Fig. 2]), 
where the stagnation point corresponds to "d c . 

4. As a consequence, no branching off takes place 
for the solutions of Bragg-Hawthorne equa- 
tion. 

5. The critical Rossby number $ c corresponding 
to stagnation is an increasing function of rj 

(Fig. ED. 

6. The whole picture can be reached by putting 
together these results with those obtained in 
PQ, where there is a branching owing to the 
boundary conditions at the frontier between 
the vortex and the irrotational flow. Moreover, 
since the results in jT] for the rotational flow 
depend on the inlet flow as well as on the ro- 
tational expansion parameter fj defined in Eq. 
(1201) . given a quasi-cylindrical approximation, 
it can be concluded that this expression is the 
minimum of v z in the core. Therefore, we can 
get the critical Rossby number $ cs and com- 
pare it with that corresponding to the fold "& c f- 
This present work verifies that $ cs < # c /, in 
accordance with Buntine and Saffman's results 
0. 

7. In the quasi-cylindrical approximation, previ- 
ous results do not depend on the chosen profile. 
This can be explained by the boundary condi- 
tions chosen depending on the inlet flow and 
on the parameter expansion. 



Acknowledgements - We would like to thank Un- 
versidad Nacional de General Sarmiento for its sup- 
port for this work, and our colleague Gabriela Di 



040002-6 



Papers in Physics, vol. 4, art. 040002 (2012) / R Gonzalez et al. 



Gesii for her advice on the English version of this 
paper. 



[1] R Gonzalez, G Sarasiia, A Costa, Kelvin waves 
with helical Beltrami flow structure, Phys. Flu- 
ids 20, 24106 (2008). 

[2] R Gonzalez, A Costa, E S Santini, On a varia- 
tional principle for Beltrami flows, Phys. Flu- 
ids 22, 74102 (2010). 

[3] R Gonzalez, E S Santini, The dynamics of 
beltramized flows and its relation with the 
Kelvin waves, J. Phys.: Conf. Scr. 296, 012024 
(2011). 

[4] The Beltrami flow is defined as a field v# 
that satisfies ug — V x vg = 7v^, with 
7 = constant. We say that the flow has a bel- 
trami flow structure when it can be put in the 
form v = Ue z + Qreg + Vg, being U and O 
constants, i.e it is the superposition of a solid 
body rotation and translation with a Beltrami 
one. For a potential flow 7 = 0. 



[5] J D Buntine, P G Saffman, Inviscid swirling 
flows and vortex breakdown, Proc. R. Soc. 
Lond. A 449, 139 (1995). 

[6] G L Brown, J M Lopez, Axisymmetric vor- 
tex breakdown Part 2. Physical mechanisms, 
J. Fluid Mech. 221, 573 (1990). 

[7] B Benjamin, Theory of the vortex breakdown 
phenomenon, J. Fluid Mech. 14, 593 (1962). 

[8] R Guarga, J Cataldo, A theoretical analysis of 
symmetry loss in high Reynolds swirling flows, 
J. Hydraulic Res. 31, 35 (1993). 

[9] S L Bragg, W R Hawthorne, Some exact solu- 
tions of the flow through annular cascade ac- 
tuator discs, J. Aero. Sci. 17, 243 (1950). 

[10] G K Batchclor, An introduction to fluids dy- 
namics, Cambridge University Press, Cam- 
bridge (1967). 

[11] S V Alcksccnko, P A Kuibin, V L Okulov, The- 
ory of concentrated vortices. An introduction, 
Springer- Verlag, Berlin Heidelberg (2007). 



040002-7 



